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'Kondo state' and Kondo resonance in a two-dimensional electron gas 



2 Stefan Miillegger,* Mohammad Rashidi, Michael Fattinger, and Reinhold Koch 

3 Institute of Semiconductor and Solid State Physics, Johannes Kepler University, Linz, Austria. 

The delicate balance of spin-screening and spin-aligning interactions determines many of the pecu- 
liar properties of dilute magnetic systems. We study a surface-supported all-organic multi-impurity 
Kondo spin system at the atomic scale by low-temperature scanning tunnelling microscopy and 
-spectroscopy. The model system consists of spin-1/2 radicals that are aligned in one-dimensional 
chains and interact via a ferromagnetic RKKY interaction mediated by the 2DEG of the supporting 
substrate. Due to the RKKY-induced enhanced depopulation of one spin-subband in the 2DEG, we 
finally succeeded to detect the so far unobserved 'Kondo state' as opposed to the well-established 
Kondo resonance. Its cloud of screening electrons, that are virtually bound to the radicals below 
the Kondo temperature, represents the extended exchange hole of the ferromagnetically polarized 
spin chain imaged here in real space. 



4 The interaction of localized electron spins of magnetic 

5 impurities with delocalized conduction electrons provides 
e a key handle for controlling the spin polarization in di- 
7 luted magnetic systems essential for future spintronics 1 ' 2 
s and molecular electronics 3 applications. Many of their 

9 unique properties are determined by the competition be- 

10 tween screening of local spins and magnetic interactions 

11 of neighbouring spins. Below a characteristic temper- 

12 ature the impurity spin may be screened by a collec- 

13 tive of conduction electrons with opposite spin - a phe- 

14 nomenon known as the Kondo effect 4 . The screening 

15 electrons form a virtual bound singlet many-body quan- 
ia turn state with the singly occupied impurity 5 denoted 
17 herein as the 'Kondo state' 3 . At higher impurity concen- 
ls trations, i. e. at sufficiently small separations, their indi- 
io vidual Kondo states may overlap, causing an additional 

20 magnetic RKKY 6 -type interaction of individual impurity 

21 spins 7 that competes with Kondo screening. 

22 Recently, the study of single impurity spins on 

23 atomic scale has become a reality by means of elec- 

24 tron transport through lateral quantum-dots 8-10 formed 

25 in a two-dimensional electron gas (2DEG) as well as 
2a surface-supported <i-metal impurities 11-13 or metalor- 

27 ganic molecules 14 ' 15 inside the tunnel junction of a scan- 

28 ning tunnelling microscope (STM). Studies of single 
20 Kondo impurities typically report a resonance observed 

30 in the electron transmission spectrum near the Fermi 

31 energy, Ep, 16 which is commonly rationalized by the 

32 Anderson single-impurity model 17 of localized magnetic 

33 states in metals [see detailed description of Fig. 1]. Al- 

34 ready in earlier work 18 the resonance signal was uti- 

35 lized for obtaining spatial information on the 'Kondo 

36 screening cloud' 19 ' 20 , i. e. the spatial region captured by 

37 the screening electrons surrounding the impurity. While 

38 the Kondo resonance involves only conduction electrons 
30 (near Ep), the virtual bound screening state (=Kondo 

40 state) is formed by all electrons screening the localized 

41 spin, i.e. the conduction electrons near Ep as well as 

42 deeper lying band electrons. A direct detection of the 

43 complete screening cloud (exchange hole) representing 

44 the Kondo state itself, has remained elusive until now. 

45 Here we report on a low-temperature STM study 

46 of surface-supported dimers and chains of the organic 




Figure 1. Kondo resonance, (a) Energy level configura- 
tion of the Anderson single- impurity model 17 . In this model 
the energy e of the magnetic impurity level (M) is assumed 
to lie sufficiently far below Ep so that its occupation proba- 
bility p is close to one, thus avoiding mixed valency (p < 1). 
Double occupation of the impurity level costs the Coulomb 
charging energy U and raises this state to M* well above Ep, 
thus making addition of a second electron energetically un- 
favourable. In the bulk of metals containing dilute magnetic 
impurities the Kondo effect gives rise to a higher density of 
states 5 for conduction electrons with opposite spin localized 
at the impurity sites. As a result the scattering cross-section 
for the conduction electrons is increased, which explains the 
resistivity increase of bulk Kondo systems below the Kondo 
temperature Tk 4 . (b) In the case of tunnelling transport the 
low-lying impurity state opens up a new transport channel 
close to Ep 3 ' 16,21 ' 22 . Whenever an electron tunnels off the 
impurity, leaving behind an empty virtual intermediate state, 
it is replaced almost instantaneously by a conduction electron 
of opposite spin as demanded by the high occupation prob- 
ability close to one. The extra transport channel increases 23 
the transmission (conductance dl /dV) close to Ep and is com- 
monly denoted Kondo resonance (R) in the literature 3,16 . 

47 sp radical BDPA (a,7-bisdiphenylene-/3-phenylallyl, 

48 C33H22) [Fig. 2a], which is a well-established spin-1/2 
40 standard in electron spin resonance spectroscopy 24 . We 
so finally succeeded to tunnel out of the Kondo state it- 

51 self and thus were able to detect and directly visual- 

52 ize the exchange hole 25,26 generated by a linear chain 

53 of molecular Kondo impurities in a 2DEG. In the present 

54 work, experimental detection of the Kondo state is en- 

55 abled by a remarkably strong change in conductance fa- 
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56 cilitated by the exceptional combination of properties 

57 of our Kondo system - as opposed to atomic d-metal 
ss Kondo systems in the literature: (i) Weak physisorp- 
59 tion of the radicals on the substrate preserves their mag- 
eo netic state, (ii) The reconstructed (111) surface of Au 
ei provides the crucial combination of surface anisotropy 

62 (herringbone reconstruction) together with a Shockley- 

63 type surface state near Ep. The former guides the self- 

64 assembly of structurally well-ordered BDPA chains and 
es the latter acts as a 2DEG capable of communicating the 
ee Kondo effect over large distances of more than a nanome- 
e? tre across the sample surface (typical for coinage metal 
es (111) surfaces 18,27 ' 28 ). Thus individual Kondo states of 

69 neighbouring BDPAs are allowed to overlap, (iii) Since 

70 the singly occupied highest molecular orbital (SOMO) 

71 of BDPA is symmetry-compatible with the sp wave func- 

72 tions of the 2DEG, constructive interaction of neighbour- 

73 ing spins is facilitated, (iv) The radical-radical separation 

74 in BDPA chains on Au(lll) is appropriate to give rise to 

75 an additional ferromagnetic 7 ' 29 RKKY-type interaction 

76 of individual BDPA spins, preferentially mediated by the 

77 2DEG, as well, (v) The resulting constructive super- 

78 position of individual Kondo states leads to an effective 

79 depopulation of one spin-subband and thus an enhance- 
so ment of the conductance change enabling detection of the 
si joint Kondo state by the STM tip. 

82 BDPA self-assembly on Au(lll) 

83 The STM topograph of Fig. 2b reveals the forma- 

84 tion of benzene-free BDPA nanostructures upon depo- 
ss sition onto Au(lll) at 300 K. Individual BDPA radi- 
se cals are imaged as bean-shaped protrusions [see also in- 
87 set of Fig. 2c] . Pronounced clustering dominates already 
as at submonolayer coverages, and dimers of BDPA consti- 

89 tute the smallest clusters. Larger BDPA nanostructures 

90 exhibit characteristic star-like shapes consisting of three 

91 branches of one-dimensional BDPA chains up to several 

92 nanometres long (depending on coverage). Single BDPAs 

93 are rarely observed on the surface due to an increased sur- 

94 face mobility even at 5 K, which strongly constrains the 

95 energy range accessible by STS compared to dimers and 
ge chains. As determined from STM topographs, BDPAs in 

97 chains pack closer (0.73 nm) than in dimers (0.8-0.9 nm). 

98 In the bulk phase BDPA is a one-dimensional paramag- 

99 net above 1.7 K forming similar linear chains with a sig- 

100 nificantly larger radical-radical spacing of 1.4 nm along 

101 the b crystallographic direction. 

102 Kondo resonance 

103 The unpaired electron of BDPA resides in the SOMO 

104 that derives from sp hybrid orbitals and is delocal- 

105 ized over the fluorene units of the molecular backbone 

106 [Fig. 2a] . When the STM tip is located over BDPA 

107 adsorbed on Au(lll), typical dl/dV spectra exhibit a 
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Figure 2. BDPA on Au(lll). (a) Chemical structure 
and DFT-calculated SOMO spin density of BDPA; the two 
'hot spots' are centred around the meso-carbons. (b) STM 
topograph of self-assembled BDPA clusters on the Au(lll) 
surface (+20 mV, 200 pA). (c) dl/dV spectrum of a single 
BDPA radical (red) and a dimer (black); the zero-bias con- 
ductance of the pristine Au(lll) background at the Fermi 
level was calibrated to about 1.3 nS; solid lines are numerical 
Fano fits of Kondo antiresonance (see also Tab. I); green ar- 
rows mark vibrational satellites; insets show STM topographs 
of single BDPA and dimer; crosses mark the STM tip posi- 
tion for recording dl /dV spectrum, (d) dl/dV Kondo spec- 
tra over different positions (o) of a BDPA dimer. (e) dl/dV 
Kondo spectra of a BDPA dimer at different temperatures, (f ) 
Temperature-dependent broadening; dashed line represents 
best-fit based on formalism of Nagaoka et al. 30 . 



TABLE I. Fano fit parameters of the BDPA monomer (aver- 
age of four), dimer (average of three), and chain (average of 
three) adsorbed on Au(lll) at 5 K; asymmetry parameter q, 
resonance position AE and half-width T. 

q (±0.1 ) AE (±0.7 meV) T (±1 meV) 

monomer —0.15 —5.9 11 

dimer -0.23 ±2.3 13 

chain -0.15 ±7.0 19 



3 



strong conductance minimum (dip) close to Ep [Fig. 2c]. 
The maximum amplitude is observed for tip positions 
near the rim of BDPA, while over intramolecular posi- 
tions the dip broadens and decreases [Fig. 2d]. The dip 
amplitude amounts to more than 20% of the zero-bias 
conductance of the pristine substrate for our best STM 
tips. We observe almost identical spectra for BDPA rad- 
icals adsorbed over different fee and elbow positions of 
the reconstructed Au(lll) surface, suggesting that the 
adsorption site relative to the substrate atomic lattice is 
negligible. A detailed analysis of our STS results confirms 
the characteristic properties of the transmission Kondo 
effect 11 ' 12 for BDPA/Au(lll): (i) The conductance min- 
imum lies within a few milli-electronvolts around E-p. (ii) 
The resonance width, 2r, increases with temperature and 
the amplitude decreases simultaneously [Figs. 2e,f[. (iii) 
The spectral shape of the dip 23 is well described by the 
Fano function 21,31 ' 32 [Fig. 2c] ; the numerical values of 
the fit parameters are listed in Tab. I. (iv) The reso- 
nance is accompanied by two satellite features symmet- 
rically offset around the Kondo dip [marked by green 
arrows in Fig. 2c] indicating molecular vibrations ex- 
cited by inelastic tunnelling processes. This is expected 
for molecule-based Kondo systems where the impurity 
orbital (SOMO) is spatially extended and thus sensi- 
tive to geometric changes of the molecular backbone 33 . 
The respective vibrational energy of rj 35-55 meV (300- 
450 cm -1 ) points to one or more collective vibrational 
modes of the BDPA backbone rather than plain C-H vi- 
brations, which are typically observed at higher energies 
as suggested by our density functional theory (DFT) cal- 
culations. 

The lower surface mobility of BDPA dimers compared 
to single radicals facilitates recording of the Kondo res- 
onance width up to temperatures of 50 K [Fig. 2f[. 
The respective broadening is analysed in detail based 
on a Fermi liquid description put forward by Nagaoka 
et al. 30 . The least -square fit based on this formalism, 
r = 2yJ (nkftT) 2 + 2(/cb2k) 2 , is shown as dashed curve 
in Fig. 2f. The obtained nominal value of Tk = 54±5 K is 
consistent with the resonance width of 2/cbIk predicted 
by Fermi liquid theory. [For further details on the Kondo 
signal of BDPA/Au(lll) see Supplementary Information 
SI.] 

The existence of a Kondo signal corroborates that 
BDPA preserves the magnetic ground state after adsorp- 
tion on Au(lll) at the single radical level. Our elec- 
tron spin resonance (ESR) experiments on samples with 
BDPA monolayer coverage on Au(lll)/mica substrates 
yield g = 1.96 at 7 K. This value is in good agreement 
with that of the single-crystal bulk phase (2.00 26) 34 as 
well as that observed in ultrathin BDPA films drop-cast 
on Au(lll) substrates (2.005) 35 at room temperature. 
Both, Kondo behaviour and ESR signal, evidence that 
a possible charge transfer between BDPA and the Au 
substrate is very small. 
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Figure 3. Spectroscopic dl /dV maps at different bias 
voltages. Overlayed yellow contour lines mark areas covered 
by BDPA molecules as extracted from STM topographs; red 
arrows indicate reduced-conductance envelopes, blue arrows 
a surface defect. 



Experimental observation of the Kondo state 

The spectroscopic images of Fig. 3 reveal intriguing 
details of the Kondo screening of BDPA clusters on 
Au(lll). Doubly occupied pure MO states [see Supple- 
mentary Information S2] are imaged as bright protrusions 
confined to the area covered by the BDPA molecules 
[Fig. 3a] . For sample bias ranging from about —0.6 to 
—0.3 V we observe areas of reduced conductance (dark) 
enveloping the BDPA clusters as marked by red arrows 
in Figs. 3c, d. These areas clearly extend more than 1 nm 
into regions of the pristine Au(lll) surface outside of 
and in between neighbouring BDPAs - in clear contrast 
to the pure MO states of Fig. 3a. We denote them as 
reduced-conductance envelopes (RCEs), since dl/dV is 
reduced by about 50% relative to zero-bias conductance 
of the pristine Au(lll) surface. The RCEs start to ap- 
pear slightly below the onset of the Au(lll) surface-state 
band close to —0.5 V [Figs. 3b, c[. They persist over an 
energy range of about 200 meV [Fig. 3d] and fade off when 
further approaching zero bias. Above —0.3 V the pristine 
standing-wave pattern of the surface state electrons is 
imaged [Fig. 3e] . It originates from scattering at the her- 
ringbone reconstruction pattern 36 of the Au(lll) surface 
and the BDPA clusters. In the standing wave pattern 
of Fig. 3e we find no evidence for an enhanced electro- 
static scattering potential localized at the BDPA radicals 
(which are hardly visible). This evidences 37 that the ad- 
sorbed radicals have (almost) zero total charge - in agree- 
ment with our ESR results. Even with the STM tip above 
BDPA the characteristic step- like conductance signal of 
the surface state survives almost unaffected [see Figure 
SI of the Supplementary Information]. This indicates 
that hybridization of BDPA orbital(s) with the surface 
state is only weak and, in particular, rules out that the 
observed RCEs originate from a possible localization 38 



4 



199 of the surface state. The possibility that local inter- 

200 ference minima of the elastically scattered surface-state 

201 electron waves would cause the observed RCEs is cer- 

202 tainly ruled out 39,40 , as well: (i) The RCEs are absent 

203 in the surface-state standing-wave pattern of Fig. 3e. (ii) 

204 The RCEs lack any oscillating behaviour expected for 

205 scattered surface-state waves 41 , (iii) The RCEs are de- 

206 tected only at the BDPA units, but not at surface im- 

207 purities [marked by a blue arrow in Figs. 3b,c,f]. Possi- 

208 ble life-time effects are known to flatten the surface-state 

209 band edge typically in the order of 10 meV 42 and thus are 

210 one order of magnitude too small to explain the observed 

211 RCEs. as well. 

212 Closely around zero bias the RCEs reappear again 

213 [Fig. 3f] and represent the conductance map of the Kondo 

214 resonance measured point-by-point in Figs. 2c ,d. The 

215 RCEs of Fig. 3f thus reflect the lateral spatial distri- 
2i8 bution of the Kondo resonance signal of BDPA clus- 

217 ters on Au(lll). The sample area at which the Kondo 

218 resonance is observed by STM has been commonly at- 

219 tributed to the spatial extent of the Kondo screening 

220 cloud 18 ' 43 . However, the interpretation of its size re- 

221 mains a controversial matter, because it is much smaller 

222 than the Kondo screening length predicted by Fermi liq- 

223 uid theory [see Supplementary Information S3]. The 

224 Kondo dip of BDPA/Au(lll) is detected more than 

225 1.5 nm laterally away from BDPA and decays according 

226 to 1 jr [see Supplementary Information S4] . This suggests 

227 that the Kondo screening of BDPA/Au(lll) strongly in- 

228 volves delocalized electronic states 22 ' 44 - most possibly 

229 the surface-state band of Au(lll) that acts as a 2DEG. 



230 Discussion: Kondo state vs. Kondo resonance 

231 The additional RCEs of Figs. 3c, d observed at ener- 

232 gies well below Ep come - at first sight - as a surprise. 

233 Since the most obvious possible alternatives are ruled 

234 out above, there is strong evidence that the additional 

235 RCEs originate from the Kondo state (KS) itself, i.e. 

236 the many-body virtual bound state including the spin- 

237 screening contributions from deeper lying band electrons 

238 as well as conduction electrons near Ep. This is illus- 

239 trated in Fig. 4: (i) Their energy onset coincides remark- 

240 ably well with the band edge of the Au(lll) surface-state. 

241 The latter is strongly involved in the Kondo screening of 

242 BDPA/Au(lll), as indicated by our results, (ii) The ge- 

243 ometric size and shape are qualitatively similar to the 

244 RCEs of the Kondo resonance, but the signal amplitude 

245 is much stronger and they extend even further into ar- 

246 eas of the pristine Au(lll) surface [compare Figs. 3d and 

247 f[. (iii) The Kondo state involves an interaction of the 

248 surface-state and the SOMO, suggesting a similar energy 

249 range for both. This is illustrated in Fig. 4 and corrob- 

250 orated by point spectroscopy [see Supplementary Infor- 

251 mation S2[. 

252 The observed reduction of the tunnelling conductance 

253 is rationalized as follows: The density-of-states of the 
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Figure 4. Schematics summarizing the dl/dV results 
of BDPA/Au(lll). In addition to the Kondo resonance 
(R) close to Ef the Kondo state (KS) is detected experimen- 
tally via contributions from deeper lying surface-state band 
electrons. The Kondo state originates from resonant spin ex- 
change processes between the magnetic impurity (SOMO of 
BDPA) and the surface-state band of Au(lll). 



2DEG consists of two overlayed identical step- functions 46 
with different Fermi wave vectors, Rashba-split in k- 
space 45 but having the same energy onset. Our results 
indicate that Kondo screening depletes one of the spin 
subbands in the surface area around BDPA (e.g. the 
up-spin subband if BDPA is up-spin). The decreased 
occupation of one subband (exchange hole) reduces the 
number of available electrons for tunnelling out of the 
surface state. The spin screening is effective not only 
for conduction electrons near Ep, but for deeper lying 
band electrons, as well. This equals an effective loss of a 
tunnelling transmission channel and thus a reduction of 
the tunnelling conductance - as evidenced by our exper- 
iments. 

The dl/dV images of Figs. 3c, d corroborate that the 
Kondo state (i. e. the screening cloud) of individual 
BDPA radicals on Au(lll) are so large that they readily 
overlap with those of their next neighbours. In the case 
of two or more neighbouring Kondo impurities an addi- 
tional RKKY-type indirect-exchange interaction of their 
spins evolves, according to predictions based on Fermi- 
liquid theory 7 . The RKKY interaction of neighbouring 
Kondo impurities is either ferromagnetic (Jrkky > 0) 
or antiferromagnetic 14 ' 29,47 ' 48 (<0) and competes with 
Kondo screening (JKondo < 0) as illustrated in Fig. 5a. 
Similar to the Kondo screening, the RKKY interaction 
between BDPA radicals in chains on Au(lll) is predom- 
inantly mediated by surface-state electrons. The sur- 
face band of Au(lll) starts 500 meV below Ep, and 
has a Fermi wave number of fcp 1 = 0.617 and 0.540 nm 
including the Rashba splitting 45 . 



A possible coupling 
mediated by bulk conduction electrons decays at much 
shorter distances due to the smaller k^ 1 = 0.08 nm 45 
and thus is negligible for intermolecular interactions in 
BDPA clusters. A calculation of the distance-dependent 
oscillation of the RKKY coupling constant for the two- 
dimensional case 49 based on the momentum and effec- 
tive mass (m* — 0.25mo) 45 of the Au(lll) surface state 
yields a ferromagnetic RKKY interaction for separations 
smaller than 1 nm. Thus, for next-neighbour BDPA rad- 
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clouds of two neighbouring BDPA radicals adsorbed on the 338 subband in the 2DEG and opens a new way to study 
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latory behaviour of the two-dimensional RKKY coupling con- 340 pled spin clusters with impact on the development of 
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for k' 1 = 0.614 nm. 341 

molecular-based magnetic logic applications for future 
342 spin- and molecular electronics. 

294 icals in chains, which are separated by 0.73 nm, a parallel 

295 spin alignment is favoured as illustrated in Fig. 5b. Ferro- 

296 magnetic RKKY facilitates a constructive superposition 

297 of individual Kondo states, which enhances the depopula- 343 Methods 

298 tion of one subband and, consequently, further lowers the 
290 conductance for tunnelling out of the surface-state. This 

300 enhancement effect seems to be essential for enabling "4 BDPA recrystallized in benzene was thermally evap- 

301 experimental detection of the Kondo state. Since the ™ orated in ultrahigh vacuum (UHV) from a quartz cru- 

302 RKKY interaction strength decays with 1/r 3 , the influ- 346 cible at 383 K after thorough degassing at 373 K. The 

303 ence of the next-neighbour BDPA is dominant in BDPA 347 single-crystal Au(lll) surface was prepared by repeated 

304 chains. Thus BDPA chains on Au(lll) are expected to "8 cycles of 0.5 keV Ar+ bombardment and annealing at 

305 be paramagnetic at 5 K in agreement with our results as "9 720 K. STM and STS experiments were carried out at 

306 well as a recent ESR study 35 . 350 7 K employing electrochemically etched W tips. The 

307 Tunneling out of the Kondo state is further supported 351 dl/dV signal was obtained from the first-harmonic cur- 

308 by the fact that the molecules are weakly imaged by STS 352 rent signal detected by lock-in technique (0.5-2 kHz; 5- 

309 between -0.7 and -0.5 V but not at -0.1 V. 353 2 mV sinusoidal peak-to-peak voltage; average of 10 

310 In dimers the radical-radical separation of 0.8-0.9 nm 3 " single spectra). W tips were deoxidized by anneal- 

311 is larger and close to the transition between ferro- and 355 in g in UHV; impurity- and tip effects were minimized 

312 antifcrromagnetic alignment [Fig. 5]. Consequently, the ™° by careful sample preparation and multiple tip form- 

313 RKKY-based enhancement effect for the Kondo state sig- » 7 ings between the dl/dV experiments, resulting in Au- 

314 nal of BDPA dimers is decreased - consistent with our re- 358 coated STM tips. Reliable tip performance was estab- 

315 suits [Fig. 3c,d]. In the case of the mobile monomers, we 359 lis hed by accurately reproducing the dl/dV signature 

316 did not succeed to detect the Kondo state, although the 380 of tne Au(lll) surface state from literature 36 . dl/dV 

317 Kondo resonance close to E F is clearly observed [Fig. 3f]. 381 spectra were recorded at constant-height conditions, the 

318 This is consistent with the absence of a RKKY-based "2 spectroscopic images (dl/dV maps) simultaneously with 

319 enhancement effect in monomers. The interaction be- 383 constant-current STM topographic imaging. Gas-phase 

320 tween neighbouring radicals in the chains is observed to 384 density functional theory (DFT) single point energy cal- 

321 affect the width of the spectral Kondo resonance. Ta- 385 culations were performed with Gaussian 03 package 50 us- 

322 ble I lists the values of the resonance width, T, for the 388 ing B3LYP hybrid functional 51 , 6-31G(d) basis set and 

323 BDPA monomer, dimer and chain as determined from 387 nxe d molecular geometry extracted from BDPA bulk 

324 our experiments. 388 crystal structure . 



Conclusions 



326 We succeeded to detect and image the joint Kondo 389 References 
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SI. Details of the Kondo resonance of 
BDPA/Au(lll) 

The shape of the Kondo resonance is determined by the 
phase between different tunneling paths (into the Fermi 
sea or into the impurity level) and the particel-hole sym- 
metry of the conduction band 1 . The observation of a dip 
(reduced conductance) rather than a peak is reflected by 
the numerical value of the asymmetry parameter q in the 
fitted Fano-lineshape. q depends on the ratio of tunnel- 
ing amplitudes for different possible tunneling paths 2,3 . 
The obtained value of q « 1 (see Table I) thus indicates 
predominant tunneling into the continuum of substrate 
states rather than the impurity state (which is here the 
SOMO) due to Coulomb repulsion. 

The energetic position, AE = E — E-p, of the Kondo 
dip is offset from the Fermi level. Both magnitude and 
sign of AE are clearly affected by the presence of a sec- 
ond neighbouring BDPA molecule in dimers [see Fig. 2c 
and Table I]. While in single BDPAs the Kondo dip lies 
below the Fermi level, it lies above it for dimers. A nega- 
tive (positive) offset, AE. indicates an occupied (empty) 
state. Fermi-liquid theory rationalizes the offset AE by 
the asymmetric level alignment of SOMO and SUMO rel- 
ative to E-p (the so-called electron-hole asymmetry 4 ) de- 
termined by the occupation number of the impurity level 
(SOMO). Without particle-hole symmetry the SOMO oc- 
cupation number becomes 7^ 1 and the Kondo peak gets 
displaced relative to the Fermi level. 4 When it is < 1 the 
level is less than half-full and the Kondo resonance lies 
above Ep and below it for > 1. 

The width of the Kondo dip increases by about 19% in 
dimers compared to singles (see Table I). 



S2. Point spectroscopy 

The energy that determines whether Kondo physics 
will be visible is /cbTk ~ 5 meV, is always smaller than 
the FWHM of the impurity level 4 (SOMO). To determine 
the energies and widths of SOMO and SUMO relative 
to the substrate Fermi level we recorded local dl/dV 
spectra of BDPA dimers with the STM tip positioned 
over BDPA near the rim (Fig. SI). The Kondo dip 
close to Ep is visible (red arrow) and significantly larger 
than the modulation ('ripples') of the surface-state band 
caused by the regular herringbone reconstruction of the 
Au(lll) 5 . Typical spectra exhibit a small filled-state res- 
onance around —0.7 to —0.9 eV and a strong empty-state 
resonance around +1.5 eV marked by black arrows. We 
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Figure SI. dl/dV spectrum with STM tip over rim of BDPA 
dimer; dark yellow: pristine Au(lll) surface. 

attribute them to the doubly occupied HOMO and dou- 
bly unoccupied LUMO-related states. 

Unfortunately, no distinct SOMO /SUMO resonances 
are observed in the spectrum. This may be due to 
the steric protection of the SOMO /SUMO, which cor- 
responds to a small tunnelling amplitude into or out of 
the SUMO/SOMO - the same reason as for the observed 
small Fano q factor (see above). Nevertheless, the posi- 
tions of SOMO/SUMO can be approximated by a com- 
parison with the spectrum of the pristine Au(lll) sur- 
face in Fig. SI. Close to about ±0.5 V the conductance 
of BDPA is reduced below that of the pristine substrate, 
suggesting two antiresonances (marked by green arrows) 
that relate to SOMO/SUMO. The respective SOMO- 
SUMO energy gap of U ~ 1 eV is consistent with the 
DFT-calculated value of U = 1.56 eV obtained for a 
single BDPA radical in the gas-phase, because a slight 
narrowing of the gap by a few tenths of electronvolts is 
commonly observed for organic molecules adsorbed on 
metal surfaces. 



S3. Kondo screening cloud 

Even the largest characteristic lengths of several 
nanometers 6,7 so far obtained from experimental data are 
much shorter than the theoretical Kondo screening length 
£k = fivp / (ksT-jg) predicted to be about 0.1-1 /im (Fermi 
velocity vp — hkp /m* and effective mass m *) 8,9 . Recent 
theoretical studies seem to be just about to resolve this 
issue, and corroborating that £k plays n o significant role 
in experimental d//dU datat, but that, rather, the ex- 
perimentally observed dl /dV features decay within a few 
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Figure S2. Distance-dependent decay of the Kondo 
signal, (a) STM topo of BDPA dimer. (b) dl/dV Kondo 
spectra recorded at different tip positions as marked by circles, 
(c) Decay of Kondo amplitude with increasing distance from 
BDPA; solid lines are fits of l/r and l/r 2 . 

Fermi wavelengths 10 . 

For the Au(lll) surface state (to* = 0.25 • to , kp = 
0.162 A" 111 ) a value of £ K = 106 nm is obtained. The 
diameter of our RCEs observed close to Ep are much 
smaller than £k- (For Au bulk electrons with vp sa 1.4 • 



10 6 m/s 12 one obtains a value of £k = 198 nm). 



S4. Spatial decay of the Kondo resonance signal 

Figure S2 illustrates the results of distance-dependent 
STS measurements along a specific high symmetry direc- 
tion of a BDPA dimer. STS spectra were recorded with 
the STM tip over positions at increasing lateral separa- 
tion from the BDPA dimer as marked by the circles in 
Fig. S2a. The respective spectra shown in Fig. S2b ex- 
hibit the characteristic Kondo dip even at large lateral 
separations of more than 1.5 nm from BDPA, while the 
amplitude of the dip decreases monotonically. The spa- 
tially decaying amplitude values are plotted in Fig. S2c 
together with the best-fit curve obtained for a decay ac- 
cording to l/r (red curve) with a reduced x 2 = 1-1. For 
comparison, we fitted a l/r 2 decay as well (blue curve), 
which results in a significantly poorer agreement with the 
experimental data (x 2 = 1.7). 
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